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ON THE ESTIMATION OF BIOLOGICAL POPULATIONS 


D. B. DeLury 
Ontario Research Foundation 


Introduction. Attempts to give mathematical formulations of bio- 
logical systems are not uncommon. However, the results of these 
studies have not found their way into biological practice to the extent 
they deserve, possibly because insufficient attention has been given to 
such questions as the experimental verification of assumptions and the 
estimation of parameters. 

It is the writer’s belief that a fusion of a mathematical develop- 
ment, based on certain assumptions, with an experimental programme 
designed to test these assumptions and to provide estimates of whatever 
parameters are involved, will often yield results of the highest impor- 
tance. Some simple examples are given in this article which discusses 
information of the type yielded by the ‘‘creel census.’’ 


Statement of the Problem. Many studies of biological populations 
involve the capture of a considerable proportion of the individuals 
in these populations. Usually the records are or could be kept in such 
a way as to show, for each day (or other period of time), the number 
of individuals caught and also the amount of ‘‘effort’’ expended in 
capturing them. Records of this kind have been used to exhibit the 
‘‘availability’’ (defined as the ratio of total catch to total effort) of 
fish in consecutive fishing seasons. It is clear, though, that a com- 
parison of availabilities cannot be translated into a comparison of 
population sizes without invoking some serious assumptions which can- 
not, as a rule, be tested. The purpose of this paper is to show that 
catch-effort records can be used to form estimates of absolute popula- 
tion sizes and to test certain important assumptions. 

It will be assumed, then, that catch and effort records are available 
for a series of consecutive time intervals. The catch for a given time 
interval, specified by ¢, will be denoted by c(#) and the corresponding 
effort by e(t). The catch per unit effort for the time interval ¢ is then 
C(t) =c(t)/e(t). Now, as a population becomes depleted, the value 
of C(t) ordinarily decreases, and the amount by which C(¢) is dimin- 
ished reflects the extent of the depletion which, in turn, depends on the 
total catch and the total effort. This qualitative statement may be 
made quantitative under suitable assumptions. This is attempted in 
the following examples. 
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Ezample 1. A Study of a Lobster Population. It may be ex- 
pected that a population of lobsters provides one of the simplest sub- 
jects for a population study, because migration is believed to be unim- 
portant, and recruitment through growth takes place only at times of 
moult. It seems not unreasonable, therefore, to make assumption (i). 

(i) The population is closed, that is, the effects of migration and 
natural mortality are negligible and the time interval is to be restricted 
so as to exclude times of moult. 

Assumption (i) alone is not strong enough to support an analysis 
of the catch-effort records, in that it supplies no information about the 
response of the population to the method of capture. Let d(t¢) rep- 
resent the proportion of the population captured during the time in- 
terval ¢t. Some assumption about the nature of d(¢) is required. 

(ii) d(t) = k(t) e(t), that is, the units of effort employed during 
interval ¢, e(¢) in number, do not compete with one another. k(#) is 
evidently the proportion of the population captured during interval ¢ 
by one unit of effort. It is proposed that k(t) be called the catchabid- 
tty and e(t) the intensity of effort. 

(iii) k(t) =k, a constant. This assumption is, of course, seriously 
open to doubt and must be tested carefully against observations before 
any conclusions are based on it. 

It is shown in Appendix I that, under these three assumptions, the 
following relations hold. 


log C(t) = log N(0)) -k E(t); 
(5) C(t) =k N(0) -k K(t), 


where E(t) and K(¢) stand for the total effort and total catch up to 
interval ¢, and N(¢) is the number of individuals in the population at 
time ¢. Now values of C(t), E(t) and K(t) can be calculated from 
the catch-effort records and plotted. If the graphs turn out to be 
acceptably straight, estimates of k and N(o) follow at once. If the 
graphs exhibit serious curvature, some modification of the assumptions 
is required. In the event that the assumptions are not contradicted by 
the data, the estimates of k and N(o) can be substituted in formula 
(2’) of Appendix I to exhibit the behaviour of the population through- 
out the whole period. 

Table I gives a day-by-day record of the catch of lobsters, in 
pounds, and the effort, measured by the number of traps fished each 

* The numbering of these equations comes from Appendix I. 


146 


, 


day. These records were made in 1944 in the Tignish area of Prince 
Edward Island. Other records, not presented here, show that the dis- 
tribution of size remains fairly constant in samples taken throughout 
the season. Thus ‘‘pounds’’ may be identified with ‘‘number of indi- 
viduals’’ without exposure to bias. 


TABLE I 
Column 1 2 3 4 5 6 7 

Date Pounds Traps C(t) x 100 K(t) MC(t)+1 E(t) 

May 2 147 200 74 0 87 0 

3 2796 3780 74 0 87 0 

4 6888 7174 96 3 98 4 

5 7723 8850 87 10 94 11 

8 5330 5793 92 18 96 20 

9 8839 9504 93 23 97 26 

10 6324 6655 95 32 98 35 

11 3569 3685 97 38 .99 42 

12 8120 8202 99 42 1.00 46 

13 8084 8585 94 50 97 54 

15 8252 9105 91 58 .96 62 

16 8411 9069 93 66 .97 72 

17 6757 7920 85 74 93 81 

18 1152 1215 95 81 98 89 

20 1500 1471 102 82 1.00 90 

22 11945 11597 103 84 1.01 91 

23 6995 8470 82 96 91 103 

24 5851 7770 75 103 88 111 

25 3221 3430 94 109 97 119 

26 6345 7970 80 113 .90 122 

27 3035 4740 83 118 92 130 

29 6271 8144 89 121 95 135 

30 5567 7965 70 128 143 

31 3017 5198 58 133 -76 151 

June 1 4559 7115 64 136 81 156 

4721 8585 55 141 74 164 

5 3613 6935 52 145 72 172 

6 473 1060 45 149 65 179 

7 928 2070 45 149 65 180 

8 2784 5725 49 150 .69 182 

9 2375 5235 45 153 65 188 

10 2640 5480 48 156 .68 193 

12 3569 8300 43 158 .63 199 

162 207 


It is convenient to define the unit of effort as 1000 traps fished for 
one day and the unit of catch as 1000 pounds. Column 4 of Table I 
shows the values of C(t), obtained by dividing the number of units 
in the daily catch by the corresponding number of units of effort. 
Column 5 contains the values of K(t), formed by accumulating the 
numbers in Column 2. Note that K(t) is the total catch up to the t® 


147 


| 
a | 
|, 
Pe 
4 
| 
{ 


day. Similarly the numbers in Column 7 give the accumulated effort 
up to the ¢™ day. 

The values of In C(t) are plotted against those of F(t) in Figure 1. 
Clearly these points cannot be considered to lie on a straight line. 
However, after May 22, the values of In C(t) decrease fairly regularly 
and it seems not unreasonable to fit a straight line to the records from 
May 23 to June 12. Therefore, for the present, the records up to May 
23 will be disregarded and totals of catch and effort calculated from 
this point on. The numbers are shown in Table II and the plotted 
points and fitted line in Figure 2. Details of the fitting are discussed 
in Appendix II. 


Ly Cle) 
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The equation of the line i is 

In C(t) = 1.78529 - 0036256 $6 

Now equation (4) is written in terms of natural logarithms. The 
equivalent form using common logarithms is 
(4’) In C(t) =In [k N(0)] —k Ime E(t). 
Hence Im [k N(o)] will be identified with 1. 18529 and kine with 
.0036256. Direct computation then yields ; 


k = .0083482 
N(o) = 112.34 units of catch = 112340 pounds 


N (0) here represents the meenend of pounds in the so coment at the 
end of May 22. . 
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Substitution of these numbers in formula (2’) provides an equation 
which describes population throughout the range May 22 to June 12. 


N(t) = 112;340 


where the E(t) values are those of Table II. In particular, putting 
E(t) = 104, the number of pounds in the population at the end of June 
12 is calculated to be 47143 pounds. Thus, according to formula (2’), 
the population decreased, from May 22 to June 12, by 65197 pounds. 
The actual catch, during the same interval, was 65964 pounds. 
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The close agreement of these two numbers reflects chiefly the good- 
ness of fit of the straight line to the plotted points and gives no direct 
support to the estimate of N(o). However, nothing in the results of 
the computations is in conflict with the assumptions according to which 
these results are obtained and interpreted. It is true that other situa- 
tions than those assumed here could lead to a linear relationship be- 
tween log C(t) and E(t), but these would require a balance among 
several variables which are unlikely to admit such balance. It is to 
be supposed, therefore, that the straightness of the line gives consider- 
able support to the assumptions on which this computation is based. 
However, no decision on the adequacy of the assumptions can be made 
with confidence until these investigations are planned with a view to 
supplying the information needed for a critical test. Good tagging 
records, for example, or a breakdown of the catches according to size 
or age, will often indicate one set of assumptions in preference to 
others. 
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TABLE II 


Date C(t) x 100 K(t) In C(t) +1 E(t) 
May 23 82 0 91 0 
24 75 7 88 8 
25 94 13 97 16 
26 80 16 90 19 
27 83 22 92 27 
29 89 25 95 32 
30 70 32 84 40 
31 58 37 -76 48 
June 1 64 40 81 53 
55 45 74 61 
5 52 50 33 69 
6 45 53 65 76 
7 45 54 65 77 
8 49 55 .69 79 
9 45 57 65 85 
10 48 60 .68 90 
12 43 62 63 96 


During the period from May 2 to May 22, the catch per unit of 
effort remains practically constant. It would seem that little can be 
gained by discussing these observations, because it is clear that catch- 
effort records in which the catch per unit effort is constant can supply 
no worth-while information about the population unless, on the basis 
of other information, this constancy can be accounted for. No such 
information is available here. 
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Instead of plotting log C(t) against E(t) the values of C(t) may 
be plotted against those of K(t). A straight line fitted to these. points 
is to be interpreted according to equation (5). When the values of 
C(t) and K(t) in Table I are plotted, a set of points much like those 
of Figure 1 is obtained. This graph is not shown here. Figure 3 
shows the values of C(t) and K(t) from Table II and the line fitted 
to them. This line has the equation 


C(t) = 0.92873 — .0079835 K (t) 
which yields the estimates 


k = 0079835, N(o) =116.33 units of catch=116330 pounds. These 
numbers agree reasonably well with those based on equation (4). 


Equation (5) is somewhat simpler for numerical work than equa- 
tion (4), because logarithms are not needed. Also, equation (5) fur- 
nishes a graphical estimate of population size, given by the intercept 
of the line on the K(t) —axis. It might be concluded, therefore, that 
equation (4) should be dropped entirely in favour of equation (5). 
However, the final decision on this matter must await the development 
of an adequate sampling theory. 

Another point which indicates that both lines should be used arises 
from the fact that, in most cases, the assumptions on which these equa- 
tions are based will not be entirely fulfilled. The inadequacy of the 
assumptions is likely to be reflected in different ways in the two equa- 
tions. It follows that substantial agreement between the two sets of 
estimates gives qualitative support to the assumptions. 

The data used in this example are not records for the whole Tig- 
nish area and therefore these population estimates do not refer to the 
whole area. It is, perhaps, worth remarking that if, in addition to the 
catch-effort records used here, a record of total daily catch or total 
daily effort were available for the whole area, then equation (5) or 
equation (4) would be fitted to obtain an estimate of the total popula- 
tion, assuming, of course, that the values of C(t) computed from the 
sample and the catchability in the sample area are representative of the 
entire area. Under slightly stronger assumptions, the total catch or 
the total effort for the whole area could be used to expand the popula- 
tion estimate for the sample area to one for the whole area. 

No evidence is available to check the accuracy of the above esti- 


mates, but they agree in a general way with what is known about this 
lobster fishery. 
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Example 2. An Application to Some Speckled Trout Records. 
Figures 4, 5 and 6 show In C(t) plotted against E(t) for the daily 
eatch-effort records made on speckled trout in Montague Pond, Prince 
Edward Island, during the years 1944, 1945 and 1946. Catch is meas- 
ured by the number of fish caught each day and the corresponding effort 
is given in rod-hours. 

The three graphs are seen to be similar in some respects. Most 
striking is the more or less steady decline of In C(t) with E(t) down 
to a well-defined minimum, followed by a pronounced upward trend. 
The day-to-day fluctuation of C(t) is very large, which seems to imply 
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that the assumption of constant catchability is not satisfied. This 
fluctuation may, however, be regarded as exhibiting non-uniformity 
in the unit of effort. In any case, the assumption of constant catch- 
ability should provide a usable description of the behaviour of the 
population provided there are no substantial trends in the values of 
the catchability. An attempt will therefore be made to account for 
the behaviour of C(t) according to the following assumptions. 


1. From the opening of the fishing season until the minimum was 


reached, recruitment was negligible and the catchability was 
constant. 


2. From the date on which the minimum was reached until the 
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end of the fishing season, the catchability was constant but the. 
population was continually being built up by recruitment. 
The selection of the date on which the minimum is reached is arbi- 
trary, but in these examples this is not serious. The dates chosen are 


shown on the graphs. 


These assumptions permit the fitting of straight lines to the first 
portions of the graphs and the estimation of N(o) and k, exactly as in 
Example 1. The results of the calculations are shown in Table III. 

The estimates for 1945 contrast sharply with those for 1944 and 
1946. The catchability in 1945 is about double those of the other 


fie. 


Serr se 
TABLE III 
Estimates Estimates Caleulated Cateh 
(equation 4) (equation 5) deerease 
N (0) k N(0) k 
1944 
May 2-Sept. 5 5761 6.874x 10-4 5801 7.256 x 10-4 3569 3471 
1945 
Apr. 17-June 24 3025 1.313 x 10-3 3266 1.214x 10-3 2183 2185 
1946 
Apr.19-Aug.6 6014 6.531 x 10-4 6287 6.207 x 104 3375 3087 


years and the population estimate is roughly one-half those of 1944 
and 1946. These discrepancies can be accounted for by the assump- 
tion that about half of the catch and the corresponding effort were not 
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reported during 1945. This question cannot be settled here, although 
probably it could be if tagging records were available, but it raises 
an interesting possibility for the use of these methods—the detection 
of poaching! Natural mortality, of course, produces the same effects 
as unreported catch. 
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Turning now to the portions of the graphs between the minimum 
and the end of the season, neither the assumptions nor the plotted 
points indicate any specific procedure as being appropriate in trying 
to account for the behaviour of €(t). The assumptions are not sharp 
enough to require that these points follow any particular curve and 
the plotted points are so widely spread out that they do not indicate 
a curve of any definite form. However, the points show clear evidence 
of trend (except in 1944, when they are too few in number to show 
anything) and since there is nothing to indicate a curvilinear relation- 
ship, a straight line has been fitted to exhibit the trend. The 1946 
records only are used in this illustration. 

The equation of the straight line fitted to the 1946 records from 
August 6 to September 14 is 


In C(t) = 0.32022 + 0.00018159 E(t). 


where E(¢) is measured from the end of August 6. Assuming that the 
eatchability has the same constant value as in the earlier part of the 
season, it follows from the development of Appendix I [formula (7) ] 
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that the population varies according to the equation 


According to this equation, the population increased from 3200 at the 
end of August 6 [equation (4) estimates this number at 2639] to 3724 
at the end of September 14. The catch caleulated from this equation 
[Appendix I, formula (8)] is 817, whereas the actual catch is 711. 
The calculated recruitment is therefore 1340. The rate of recruitment 
is calculated to be about 1.2 per cent. 

This treatment of the second portion of the fishing season is most 
unsatisfactory, based as it is on a set of speculations which cannot be 
tested with the available data, and is given here only to show one of 
the several types of computation which might be employed in such 
cases. 

These data do not, in themselves, give any assurance that the up- 
ward trend in the log C(t) — E(t) curve is the result of recruitment 
rather than, for example, an increase in the ecatchability, although 
general experience does support this assumption. This point could be 
settled if the population had contained enough tagged individuals to 
show up well in the catches at the time when the upward trend began. 
Granted that recruitment does take place, it may be important to know 
if it is supplied by immigration or through growth. Presumably some 
information on this question would be obtained by breaking down the 
catches according to size. Indeed, with these two additions to the 
catch-effort data, many sources of perplexity would be removed. 

The one possibility of an independent check on these population 
estimates lies in some tagging records which were taken along with the 
1944 catch-effort data. On May 15, 1944, 178 marked trout were re- 
leased in the pond and all recaptures were recorded. The recapture 
records for each day have been used to estimate the size of the popu- 
lation on that day. The results of these calculations are exhibited in 
Figure 7. The solid line represents the size of the population com- 
puted from the population equation. 

The population estimates based on the recapture of marked fish have 
very large sampling errors, since each one is based on a small number 
of recaptures. This is reflected in the large dispersion of the plotted 
points, but even so, certain facts emerge clearly. The population esti- 
mates based on recaptures of marked fish show a definite downward 
trend as the season progresses and these estimates are persistently 
lower than those given by the population equation. 


155 


a 
— 
| 
| an 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
; 
| 
| 
| 
| 
| 
1 
| 
| 
| 
| 
: 
| 
| 
| 
| 
| 
| 
| 


There is no reason to regard the estimates obtained from the marked 
fish as more trustworthy than those calculated from the population 
equation. Indeed, most of them are absurdly low, in view of the num- 
ber of fish taken from the pond during this period. In any event, it is 
desirable to account for the discrepancy between the estimates given by 
the two methods. 

The bias usually feared in tagging experiments arises from the pos- 
sibility that a number of the marked individuals die as a result of the 
marking. It is not known whether or not any of these marked fish 


} 


died, but in any case, this would result in estimates that are too large. 
A bias of another sort is needed to account for estimates that are per- 
sistently low. 

If, for any reason, the marked fish are more likely to be caught 
than are the unmarked fish, population estimates based on their recap- 
ture will be too low. The assumption of such an effect may seem to 
be unwarranted, but in the case under consideration, something of this 
kind may well have happened. The fish to be marked were caught 
by the same method as was used to capture them later on and if some 
of the fish were more ‘‘catchable’’ than others, the group of marked 
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fish presumably contained a larger proportion of ‘‘highly catchable’’ 
fish than did the population as a whole. 

With this possibility in mind, it is of interest to apply to the sub- 
population of 178 marked fish the same methods of estimation as were 
used for the whole population, even though it is not to be expected that 
the methods would perform well because of the smallness of all the 
numbers involved. On a number of days, no marked fish were caught, 
which makes it necessary to group together sets of consecutive records 
so as to avoid zero values of C(t). 

The following estimates are obtained by fitting equations (4) and 
(5) to these grouped data. 


Equation (4) : (0) = 137, k = 1.0406 x 10 
Equation (5) : N(o) = 145, k = 1.1341 x 10° 


It is particularly interesting to find that the catchability of marked 
fish is almost double that of the population as a whole. Assuming that 
this is not accidental, it may well be simply a size effect and might 
disappear if these computations were performed for fish of a given size. 


Summary and Conclusions. An outline is given of a procedure for 
obtaining, from catch-effort records, estimates of the population size 
and certain important parameters. The method depends on the selec- 
tion of appropriate assumptions to account for the observed behaviour 
of the population. The choice of assumptions may often be difficult to 
make, inasmuch as several different assumptions may provide reason- 
able explanations of the data. It has been suggested in this article 
that population investigations should be so designed as to furnish infor- 
mation which will distinguish among the effects of the several plausible 
assumptions. 

Probably biologists will contend that the assumptions made in the 
examples of this article are too simple and cannot hold even approxi- 
mately in the populations which have been discussed. This may well 
be the case. For example, catchability undoubtedly varies with size 
and, in some instances, with other factors as well. Thus, the constant 
eatchability used in these examples must be regarded as an average 
over all sizes which appear in the catch. This average will be useful 
only if the distribution of size in the population remains constant, a 
condition that usually will not be met. This particular difficulty can 
be overcome by applying the methods to restricted size ranges, which 
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requires that the catches be recorded according to size. This device 
might well lead to a picture of the distribution of size in the popula- 
tion that is fairly free from bias, which cannot be accomplished (in 
the case of fishes, at least) by direct sampling. It is quite possible, 
therefore, that the proper use of stratification in an investigation will 
produce data to which simple assumptions may safely be applied. 

The assumption of constant catchability is not essential to the appli- 
cation of these methods, although, of course, the treatment is simpler 
when it can be made. Methods of this kind can be extended and 
elaborated indefinitely, but are likely to have little value unless the 
designs of the experimental investigations are correspondingly refined. 
It is obvious, too, that the successful application of these methods re- 
quires a thorough understanding of the biological factors which in- 
fluence the behaviour of the subjects of the study. 

It should be worth while, during the course of a population inves- 
tigation, to plot C(t) against K(t) and log C(t) against E(t) as the 
records are taken, because these points will give immediate indication 
of a change in the behaviour of the population and steps can be taken 
at once to investigate the causes and character of the change. 
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APPENDIX I 


Notation and a General Equation. 

N(t) represents the number of individuals in the population at 
time t¢. 

r(t) denotes the relative rate at which individuals are added to 
the population at time ¢. Thus, the number of individuals 
added to the population during the interval (t, t+dt) is 
given to the first order by 

N(t) r(t) dt. 
The function r(t) is intended to cover such sources of varia- 
tion as migration, growth, ‘‘natural’’ mortality, ete. and may 
take both positive and negative values. 

d(t) represents the relative rate at which individuals are removed 
from the population by the sampling method. 
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In terms of these symbols, the population increases, during the 
interval (t, t+ dt), by 
N(t) [r(t) -—d(t)] dt 
individuals. This may be equated to 
N(t+dt)-N(t). 


Dividing by dt and passing to the limit as dt approaches zero, this 
equality yields the relation 


[r(t)-a(t)], or 
(1) © N(t) -r(t) -d(t). 


All the equations used in this paper are derived from special cases of 
equation (1). 


Some Special Cases. The developments of example 1 and part of ex- 
ample 2 rest on the following assumptions. 


(i) r(t) =0; 


(ii) d(t) = k(t) e(t), where e(t) represents the rate of expendi- 
ture of effort at time t; 


(iii) k(t) =k, a constant. 


Equation (1) becomes, under these three assumptions, 


log N(t) =—ke(t), 


which can be integrated to give 
(2) log N(t)-log f'e(t) dt =-k E(t), 


where E(t) stands for the total effort expended during the inter- 
val (0, t). The time origin may, of course, be put in any con- 
venient place. Equation (2) may be written more compactly in 


the form 
(2’) N(t) =—kN(o)e* E(t) 
Now, from either of the equations (1’) and (2’), it follows that 
(3) GN __ kN (o)e-* 


dE 
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(4) 


(5) 


N 
The derivative ae which represents the rate of increase of the 


population per unit of effort expended, can be identified with the 
negative of the catch per unit effort, in view of assumption (i). 
Using C(t) to denote the catch per unit effort and writing — C(t) 


for ll equation (3) becomes on taking logarithms, 
dE’ “4 


log C(t) = log [k N(o)] E(t). 


Under assumption (i), the total catch during the interval 
(0, t), denoted K(t), is given by 


N(o0) -N(t) = N(o) —N(o0) e* , by equation (2’) 
Also, equation (3) under assumption (i) states that 
C(t) =kN(0) 


Elimination of E(t) between these two equations yields the rela- 
tion 


C(t) =k N(o0)-kK(t). 


Thus, under assumptions (i), (ii), and (iii), equation (4) holds and 


implies equation (5). Likewise, under the same assumptions, equation 
(5) is easily shown to imply equation (4). Therefore, provided the 
three assumptions are strictly fulfilled, equations (4) and (5) are 
wholly equivalent. 


{6) 


(7) 


It may be noted here that the relation 

C(t) =k(t) N(t) 
holds under only the rather mild assumption (ii). Equation (5) 
is simply one form of equation (6) when assumptions (i) and (iii) 
also apply. 

When assumption (i) is dropped, as it was in the attempt to 
account for the behaviour of C(t) during the later parts of the 
seasons in example 2, equations (1’) to (5) no longer apply. If 
a linear equation, log C(t) =a+b E(t) say, does describe ade- 
quately the relation between log C(t) and E(t), it follows from 
equation (6) that 


k(t) N(t) =e? 
In equation (7), values of a and b are determined by the data, 


but the function k(t) is unknown. A value was assigned to k(t) 
in example 2 by assuming, not only that k(t) was constant, but 
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that its value was the same during the later part of the season 
as it had been during the earlier part. No particular justification 
can be offered for these assumptions, but this is the sort of ques- 
tion that can be provided for in the design of an investigation. 

After a constant value has been assigned to k(t), equation (7) 
furnishes an explicit formula for N(t) : 


(7’) N(t) = E(t) 


If this formula gives an adequate description of the population, 
several computations are of interest. 

The total catch during the period in which (7’) applies is 
given by the following integral, evaluated over the whole of this 


period. 

(8) fern e(t) dt = Bt), 
where E(t) here represents the total effort applied during this 
period. 


The rate of recruitment can be calculated as follows. Equa- 
tion (1) becomes, under the assumptions in use here, 


d 


(9) Nt) = r(t) -ke(t), 
and from equation (7’), 
(10) log N(t) +log k=a+b E(t). 


Differentiating with respect to ¢t in equation (10) and using this 


result to eliminate £ log N(t) from (9), it is found that 


(11) r(t) = (b+k) e(t). 


It is of course, unreasonable to express the rate of recruitment as 
a function of the intensity of fishing. However, e(t), in spité of 
fluctuations from day to day, contributes to the equations in a 
virtually constant manner, and therefore r(t) is effectively con- 


stant. 
The total recruitment for this period is given by 
(12) fro N(t)dt = (b +k) fro e(t)dt 
= (1 4 4) times the total catch. 
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Clearly these computations can have little value in the present ex- 
amples, because they are based upon inadequate information. It is to 
be expected, however, that when the requisite information is available, 
calculations of this nature may have some importance. 


APPENDIX II 


The Fitting of the Straight Lines in Examples 1 and 2. The symbol 
C(t), as used in Appendix I, represents an instantaneous rate, the rate 
of capture of each unit of effort at the instant ¢t. Values of C(t) have 
been approximated from the data by means of ratios of differences, 
that is, (daily catch) /(daily effort), which may be regarded as daily 
averages of the values of C(t). The values of K(t) and E(t) which 
correspond to these average values of C(t) should presumably be caleu- 
lated up to points located somewhere between the beginnings and the 
ends of the days, but in practice, this is impossible. The catch and the 
effort corresponding to the value of C(t) for a given day must be 
accumulated up to that day, or up to and including that day. The 
first of these alternatives has been used in this paper. 

The choice of an appropriate method of fitting equations (4) and 
(5) cannot be made satisfactorily until a statistical model is con- 
structed to give an adequate description of the behaviour of the popu- 
lation. Lacking such a model, equations (4) and (5) have been fitted 
as if they were regression equations, with E(t) and K(t) treated as 
independent variables and In C(t) and' C(t) as dependent variables. 
This procedure appears to be not unreasonable, except that the question 
of weighting deserves attention. Some discussion of this point is given 
in Appendix III. 

APPENDIX III 


A Simple Statistical Model. A population of N individuals may be 
represented by a set of N identical white beads in a box. A unit of 
effort may be introduced as a dip with a small scoop, which removes a 
number of the beads. These beads are counted and replaced by an 
equal number of red beads. The red and white beads are thoroughly 
mixed in the box and another scoopful of beads is removed. Any red 
beads taken by the scoop are returned to the box; the white beads are 
counted and replaced by an equal number of red beads. This process 
is continued. 

The red beads are not considered to belong to the population, but 
are used only as ‘‘fillers’’, to ensure that a unit of effort always takes 
the same proportion of the population. The red beads could, of course, 
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be regarded as tagged members of the population and population esti- 
mates could be based on the proportions of red beads in the samples. 

Two sources of ‘‘error’’ may be expected in this sampling scheme. 

(i) The proportion of white beads in the sample will usually be dif- 
ferent from the proportion of white beads in the population. Dis- 
crepancies arising from this source will be called ‘‘sampling errors’’. 

(ii) The total number of red and white beads taken will vary from 
one sample to another. This may be regarded either as a variation in 
eatchability or as a variation in the unit of effort. Indeed, these two 
statements are identical. The two sources of variation (i) and (ii) 
will be considered separately, even though they are not wholly inde- 
pendent. 


(i) Assume for the present that each dip with the scoop takes out 
red and white beads, so that the catchability is strictly constant. Let 
the proportion of white beads in the box just before the t® sample is 
taken be p;. Assume further that the ratio of » to N is so small that 
the distribution of white beads in samples of n may be closely approxi- 
mated by a binomial distribution. Then, if C(t) is the number of 
white beads taken in the t” sample, the mean value and variance of 
C(t) are 


&C(t) =n pt, Var C(t) =n 


If K(t) white beads have been removed in the first (t-—1) samples, 
then p; = 1— K(t)/N and the mean value of ((t), conditional upon the 
total number of white beads taken in the first (f -1) sampies, is 


Evidently, n/N is the catchability, k, and therefore 
(5’) & (C(t) | K(t)] = kN-k K(t). 


This is the form taken by equation (5) for this statistical model. 

Since no question of the sampling behaviour of K(t) enters into 
equation (5’), it follows from least squares theory that the linear esti- 
mates of the parameters (kN) and k which have minimum variance 
are those which minimize 


Sow [C(t) -kN+kK(t)]?, 
t 


where the weighting factors w(t) are inversely proportional to 
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Thus the estimation of (kN) and k involves only the fitting of a 
weighted regression of C(t) on K(t). 

The relation of C(t) to E(t), analogous to equation (4), seems to 
be less obvious. It is clear, though, that such a relation must have a 
somewhat different character from that of C(t) to K(t), because the 
sampling behaviour of C(t), given K(t), depends only on the value 
of K(t), whereas that of C(t), given E(t), must involve the whole 
stochastic process up to the t sample. For this reason, it is necessary 
to compute the unconditional mean of C(t). The following relations 
yield this unconditional mean. 

The conditional distribution of C(t) given K(t) is specified by 


P[C(t) =2| K(t)] = Ea 


where (") = Therefore, 


P[C(t+1) =2|K(t)]= 1) =2, C(t) =a| K(t)] 


A=0 


Using this relation, 


[C(¢+1) | K(t)] = x P[C(t+1) K(t)] 


z=0 


is calculated, by direct evaluation of this sum, to be 
n(1-5) [1-22 |- KH] G-HE [CCH | 


A continuation of this type of calculation yields the relation 


A(C(t+y) | K(t)]=(1-k)Y (C(t) | K(#)], 
which holds for ¢ = 1, 2, 3,....and y =0, 1, 2, 3,.... 


When ¢ =1, K(t) =0, & [C(1) | K(1)] =n=kN 
and & (C(y+1) | K(1)] =kn(1-k)* 
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is the unconditional mean of C(y+1). Therefore, 
GC (t) = kN(1-k)* = kN(1-k)5™, 
since E(t) =t-1 in this case. Taking logarithms, 
log & C(t) =log (kN) + E(t) log (1-4). 


This is the equivalent of equation (4). log (1-—k) may be replaced 
by —k without sensible error, since k has been assumed to be very small. 
The fitting of this equation to observations falls outside the scope 
of regression theory. The line of best fit, in the least squares sense, 
is obtained by the same calculation as would be used in fitting the 
weighted regression of log C(t) on E(t), since E(t) is known without 
error, but questions involving the estimation of parameters are com- 
plicated. The weights to be used are inversely proportional to 


K(t)/[N-K(t)] 


(ii) When the catchability varies from one sample to another, as 
would be the case in this model when the total number of beads taken 
by the scoop varies from sample to sample, two cases may be considered. 

In sampling from the population of beads, the value of n, the total 
number of beads, could be recorded for each sample and exact allow- 
ance could be made for variations in the values of n. No new question 
is involved in this case and nothing of interest can be found, because 
in practice nothing equivalent to the value of can be observed. 

When values of n or their equivalents cannot be recorded, it seems 
not unreasonable to treat as a random variable, with mean n, and 
variance o,”._ Then the value of C(t), given K(t), can be expressed in 
the form 


C(t) = (1- ay) 


where 8; and A; are random variables with zero means. These variables 
are not independent, but their covariance is zero. Their variances are 


on? and n, € - Hence, 


Var (C(t) | K(t)] = +m, 20) 
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The weighting factors, therefore, are proportional to 


(13) w(t) = A= M/on? 


The value of o,? or its equivalent would not ordinarily be known 
and the values of A could conceivably range from —1 to +0. When 
A =-1, the variation in catchability is so large as to completely over- 
shadow the sampling errors; A =+ 0 corresponds to case (i) discussed 
above and would not occur in practice. It may be guessed that, as a 
rule, A will lie between -1 and +1. A takes the value +1 if n follows 
a symmetrical binomial distribution and vanishes if n follows a Poisson 
distribution. Furthermore, the dependence of the weighting factors on 
the value of A becomes important only after the ratio K(t)/N reaches 
a sizable magnitude. Therefore in many, perhaps most, cases, weights 
proportional to 


(14) (1- 


will serve reasonably well. 
The corresponding weights for fitting log C(t) are given by 


The value of N to use in the weighting formula would not be known 
in practice, but its value could be approximated by an unweighted fit- 
ting or graphically by a line fitted freehand to the (C(t), K(t)) points. 
Further accuracy could be reached by using the value of N, obtained 
from the weighted fitting, to compute adjusted weights, which would 
then be used for a second weighted fitting. This iterative process can 
be continued, but such attempts at precision would not generally be 
required. Indeed, estimates based on an unweighted fitting should be 
good enough for many purposes. 

For this statistical model, estimates of N and k based on equation 
(5’) can be supplied with standard errors by ordinary regression 
theory. If the fitted equation is written C(t) =a+6 K(t) and if ¢,,, 
C12, Coe denote the elements of the matrix inverse to that of the normal 
equations which determine a and b, then the estimated variances of 
a and b are given by s*c,, and s*c.., where n s? is the residual sum of 
squares and +2 is the number of pairs of observations to which the 
equation is fitted. Their estimated covariance is s¢,2. 
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The estimates of k and N areb anda/b. The standard error of the 


estimate of k is therefore s\/c.. and a first order approximation to that 
of N is given by 


b 


If normal theory is assumed to apply here, fiducial limits may be 
attached to k in the ordinary way, that is, by calculating 


V C22, 


where ¢ is taken at the desired confidence level with m degrees of free- 
dom. Fiducial limits for N are given by the roots of the quadratic 
equation 


N?(b? — + 2N (ab — t?s?c,2) + (a? t#s?c,,) =0 


provided b? — t?s?c.. > 0. 

Even though the distribution of C(t) may depart somewhat from 
normality, these fiducial limits should be useful. 

In sampling from this statistical model, each sample is based on 
one unit of effort, a restriction which is likely to be seriously violated 
in practice. The disturbance thus introduced into the weights could 
be removed, either by grouping together sets of consecutive records into 
groups containing approximately equal amounts of effort, or by re- 
placing the weighting factors w(t) by e(t) w(t). 

The conclusions reached for this statistical model may be summed 
up as follows: estimates of N and k should be calculated from equation 
(5), fitted using suitable weights, rather than from equation (4). 
Standard errors of these estimates may then be calculated. 

Doubtless this model ignores many factors which are important 
in natural populations. Much work, both theoretical and experi- 
mental, is needed to provide acceptable models, not only for the meth- 
ods suggested here, but for all investigations into biological populations. 
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PROCEEDINGS OF THE FIRST INTERNATIONAL 
BIOMETRIC CONFERENCE 


Marine Biological Laboratory, Woods Hole, Mass. 
Friday Morning, September 5, 1947 


As Chairman of the Organizing Committee, Dr. C. I. Bliss opened 
the meeting and called on Dr. Charles Packard, Director of the Marine 
Biological Laboratory, who spoke as follows: 


WELCOME TO THE MARINE BIOLOGICAL LABORATORY 
CHARLEs PackarD, Director 


It is a pleasure to welcome, on behalf of the Marine Biological Labo- 
ratory, the First International Biometric Conference and to greet old 
friends and new who are attending it. We are glad to have meetings 
of this kind here in Woods Hole for they bring scientists whose primary 
interest in the general field of Biology are somewhat different from our 
own. It may be that your presence here will stimulate our investiga- 
tors to apply those biometric principles which will aid them in formu- 
lating and interpreting their results. 

When this Laboratory was founded 60 years ago the Director, Dr. 
Whitman, looked forward to the time when all kinds of biological re- 
search should be carried on here; not only morphology and embryology 
which were then the chief topics, but also less developed fields including 
Bacteriology and Plant and Animal Breeding. Were he alive today 
he would certainly put Biometrics high in the list of desirable studies, 
for much of the work here now is quantitative in nature. 

This Laboratory has grown from small beginnings in a little 
wooden building where a few investigators and students gathered each 
summer, until now it can accommodate 500 investigators, assistants, 
and students. With this growth has ceme specialization. We have a 
great library, probably the most useful ‘and the most used library to be 
found anywhere. Our Apparatus and Chemical Departments provide 
our workers with all the usual tools of research; the Supply Depart- 
ment brings in quantities of material every day for the investigators 
and students; we have also a machine shop, a glass blowing shop, and a 
microfilm service. All of these departments are open for your in- 
spection. 

This is a democratic community where students and young inves- 
tigators feel free to ask advice from more experienced workers, and to 
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meet them informally at the Mess and on the bathing beach. This is 
one of the delightful aspects of life here, and is a reason why our people 
return year after year. We hope that the members of this Conference 
also will frequently come back to Woods Hole. 


Dr. Bliss then reported that he had asked Doctors A. Linder, Carl 
F. Kossack and E. W. Lindstrom to serve as a committee to nominate 
permanent officers of the Conference. This committee proposed the 
following nominees: Chairman, Georges Teissier; Co-Chairman, R. A. 
Fisher ; Secretary, C.J. Bliss. Dr. N. Rashevsky seconded the nomina- 
tions and the above officers were elected. 

In the absence of Dr. Teissier, Professor Fisher took the chair and 
called upon Dr. Bliss for the report of the Organizing Committee. 


REPORT OF THE ORGANIZING COMMITTEE 


The first step leading to the present Conference was the appoint- 
ment on April 8 of a committee of four by D. B. DeLury, Chairman of 
the Biometrics Section of the American Statistical Association. Dr. 
DeLury noted that the International Statistical Conferences in Wash- 
ington made no specific provision for biometry. Since statisticians 
interested in biological statistics would be in America for the Wash- 
ington meetings, he asked the committee to arrange a conference at 
which international cooperation in biometry could be discussed and 
placed on an effective basis. In accepting this assignment, the Organ- 
izing Committee added to its number so as to be broadly representative 
of quantitative biology in its statistical and mathematical aspects. 
Two representatives were named by the National Research Council. 

Financial support was provided by the Rockefeller Foundation 
which contributed $1,000 to meet the travelling expenses of delegates 
from overseas between New York and Woods Hole and other costs. 
The Marine Biological Laboratory agreed. to receive these funds and 
disburse them on order of the Organizing Committee. Any surplus 
reverts to the Rockefeller Foundation. , 

In preparing the list of foreign invitees to the Conference, we had 
the active cooperation of the Joint Arrangements Committee for the 
International Statistical Conferences in Washington. Invitations were 
sent to those planning to attend the sessions in Washington, who we 
believed might be interested in this Conference. To this list we added 
the names of other prominent foreign biometricians, who were then 
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invited to the Washington meetings by the Joint Arrangements Com- 
mittee. Invitations to the Biometric Conference were sent to 52 scien- 
tists representing 18 different countries. The Organizing Committee 
also sent individual letters of invitation to 157 American and Cana- 
dian scientists and announced the Conference in Science and elsewhere. 

A major objective of the present Conference is to lay the founda- 
tion for future cooperation in biometry. International cooperation 
in other sciences has taken various forms, ranging from triennial con- 
gresses arranged by small committees appointed by the preceding con- 
gress to formal organizations with individual scientists as members. 
Dr. Stuart Mudd, of the International Union of Biological Sciences, 
has urged the advantages of formal organization and recognition by 
UNESCO, which would make us eligible for financial aid in respect to 
future meetings. A letter from Dr. Joseph Needham, Head of the 
Division of Natural Sciences in UNESCO, expressed regret that he 
could not attend personally, wished the Conference every success and 
hoped it would be the first of many. He has also asked whether we 
‘‘intend to set up an international committee or society of a perma- 
nent nature to concern itself with ihis subject (biometry) and these 
conferences.’’ A letter from Dr. Stuart Rice notes that the proposed 
reorganization of the International Statistical Institute provides for 
affiliation by international organizations in specialized fields of appli- 
eation. Saturday morning has been set aside for discussing this prob- 
lem and for forming some continuing international biometric organi- 
zation if that is the wish of the Conference. 

The Organizing Committee is especially indebted to our host, the 
Marine Biological Laboratory. The arrangements which have been 
made with its cooperation are indicated on the printed program, to- 
gether with a brief sketch of the biological institutions in Woods Hole 
which make it a major center of biological activity in the United States, 
and a logical choice for the present Conference. 

Respectfully submitted, 


The Organizing Committee 
C. I. Bliss, (Chairman), Connec- E. J. deBeer, The Wellcome Re- 


ticut Agricultural Experiment search Laboratories 
Station H. K. Hartline, University of 
A. F. Blakeslee, Smith College Pennsylvania School of Medi- 
W. G. Cochran, North Carolina cine 
State College H. W. Norton, Oak Ridge, Tenn. 
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John von Neumann, Institute for E. W. Sinnott, Yale University 


Advanced Study J. W. Tukey, Princeton Univer- 
N. Rashevsky, University of Chi- sity 
cago E. B. Wilson, Harvard University 


The report was accepted. Dr. Bliss moved that the Chairman name a 
committee to consider what type of international organization would 
be most suitable for biometry at the present time, to submit for con- 
sideration a draft constitution embodying its recommendations and 
to report to the Conference at the morning session on September 6. 
The motion was carried. 

The meeting was then turned over to the chairman of the scientific 
session, Dr. A. F. Blakeslee, who called on Professor Fisher for his 
paper on ‘‘A Quantitative Theory of Genetic Recombination.’’*® 


Friday afternoon, September 5 


The Chairman, R. A. Fisher, proposed the following committee 
to consider the form of international cooperation: Doctors Belz, Bliss, 
Bose, Bronk, Dieulefait, Fisher, Hopkins, Linder, Neurdenburg, 
Rasch, Teissier and Tukey. This proposal was adopted. A Commit- 
tee on Resolutions was named consisting of D. G. Catcheside, Jacques 
Monod and E. B. Wilson. 

Professor Wilson then took the chair and called on several speak- 
ers to report on recent biometric developments overseas. 


* This paper, and the ensuing discussion, will be published in the March 1948 issue 
of Biometrics. 
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RECENT BIOMETRIC DEVELOPMENTS IN DENMARK 
G. Rascu 


In reporting informally on recent developments of biometry in 
Denmark, I shall review briefly some of the problems which we have 
studied. Because they are more familiar to me, I shall deal especially 
with the investigations of my colleagues and myself. The extent to 
which our results are new or merely corroborative of similar studies 
in other countries is a question I shall leave to others. The first three 
sections deal primarily with biomathematical problems and the re- 
mainder with statistical aspects of biometry. 


BIOMATHEMATICS 


Descriptive or empirical relations. Th. Madsen, partly in collabo- 
ration with S. Arrhenius, has described several immunological proc- 
esses in terms resembling mono- and bimolecular chemical reactions. 
This work was continued by Jensen [5] and in 1941 Ipsen [2] pub- 
lished a comprehensive study of the haemolysis of red blood cells and 
the length of survival of mice, rabbits and other experimental animals 
following different concentrations of several toxins. He was able to 
reduce both types of response to an equation of the form 


where D is the dose producing a 50 per cent response, 7 the corre- 
sponding reaction time, d and ¢ are the asymptotes for dose and time, 
a and K are constants. 

In the field of growth the publication of Wetzel’s control chart for 
the growth of children [13] led me to analyze the growth in height and 
weight of 150 Danish school children. When plotted against each 
other, the logarithms of height and weight for each individual at suc- 
cessive age intervals followed a straight line very closely. The lines 
varied considerably from one child to another, both in position and in 
slope. The slopes in particular varied much more than was appar- 
ently the case for the American children on whom Wetzel based his 
growth channels [11]. 

In analyzing individual growth curves more generally, 1 have 
found a relation which was applicable to data on bacteria, mice, rats, 
calves and babies. The logarithm of weight can be plotted as a 
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straight line against the age measured in a metric characteristic of each 
species. This metric or ‘‘ physiological age’’ is determined empirically. 
While varying considerably between species, it is applicable to normal 
individuals within a species {12}. 

Deductive. An interesting theory on contact epidemics has been 
developed by Petersen [9]. From general assumptions on the fre- 
quency of chance contacts of an infected with an uninfected person, he 
derived a differential equation for the spreading of an epidemic. On 
integration this led to the differentiated logistic curve, which of course 
could be discriminated from other similar curves only in very large 
epidemics. But his underlying assumptions specified the dependence 
of the duration of an epidemic on the number of cases and on the size 
of the population. From the data of various epidemics (scarlatina, in- 
fluenza, etc.) in towns differing in size from 3000 to 600,000 inhabitants, 
he found this relation to be satisfied. 

Other deductive relations, applying to laboratory studies, must be 
omitted here for lack of time. 

Biological standardization. Some of my first experiences in bio- 
logical assays drew my attention to an inconsistency between what bio- 
logical standardization professed to do and the results of experiments. 

The fact had emerged from early experimental work that absolute 
measurements of a drug in ‘‘frog units’’ ete. were impossible due to the 
varying response of different strains of the same species of animal. 
The introduction of standards was intended to remove this source of 
error by making each measurement a comparison with a standard - 
preparation, preferably an international one. To each unknown 
preparation was ascribed a ‘‘relative potency’’ expressed in terms of 
(international) units of the standard. 

It follows from well known considerations that if an wnambiguous 
relative potency can be derived from an experiment, then the reaction 
curves as plotted against the logarithms of the dose must be parallel. 
Experimental evidence is now available from widely different fields 
which shows that this condition is not always fulfilled (various sera, 
thyroxins, A-vitamins, tuberculins, ete.). 

The problem, however, goes much deeper. Even when the reaction 
curves are parallel, it is clear that an unambiguous statement about the 
relative potency of two substances requires that we get the same answer 
from different species of experimental animals as well as from differ- 
ent techniques applied to the same species. Neither of these condi- 
tions holds generally. Despite parallel reaction curves, Pedersen- 
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Bjergaard [8] has shown that the ‘‘relative potency’’ of estrone and 
of estradiol, for example, as determined in mice, rats and guinea-pigs 
varied from about 1:10 to 10:1. Similar results are known for peni- 
cillin. 

The chemistry of antibiotics, hormones and vitamins is now so far 
advanced that the explanation of these facts is obvious. The prepara- 
tions compared consist of several chemical compounds, all biologically 
active but in different ways. If the compounds in two preparations 
were present in unequal proportions, all sorts of results could come out 
of various kinds of comparisons. 

The standardization of several medicaments is now chemical, but 
for many others we still have to rely on biological methods. The pres- 
ent criticism does not preclude the possibility of a biological control 
but I believe that the methods of biological standardization should be 
reconsidered. A single figure can no longer be taken as the relative 
potency of a substance unless it is chemically identical with the 
standard. 


BIOSTATISTICS 


Distributions. We have found it useful at the State Serum Insti- 
tute to study the distribution of biological data before applying the 
usual tests. One reason is technical. Consider, for example, a skew 
distribution which on some transformation (to logarithms, square roots 
or anything else) is normal. The specific parameters of the trans- 
formed distribution are the mean and the variance, but the mean and 
the variance of the original distribution are rather complicated func- 
tions of these specific parameters. Hence, if two sets of original obser- 
vations are compared by the usual methods based upon the normal dis- 
tribution, characteristics may be lost which are obvious after trans- 
formation. 

The distribution as a whole is considered to be an important bio- 
logical datum in itself, which, in time, should be derived from general 
principles. Kapteyn’s distributions [6] have been derived from a gen- 
eral, although rather simple principle, viz., from stochastic differential 
equations. These are supposed to be approximate descriptions of the 
stochastic processes generating the distributions. We have been col- 
lecting data of different kinds which on transformation give normal 
distributions, with the intention of making a comprehensive study of 
these transformations. We would be especially grateful for sugges- 
tions which would enable us to enlarge this collection. 
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Subnormal variation. In microbiology it sometimes happens—as 
shown by Berkson in tlhe case of erythrocyte counts—that the cells 
within a preparation are distributed more evenly than would be ex- 
pected if they were distributed independently and at random. A strik- 
ing example of reticulocyte counts was reported in Denmark [10]. 
Using a meticulous technique, C. M. Plum obtained routine counts from 
different suspensions of the same blood sample. On statistical analysis 
these showed a variance of about 9, independent of the mean count. 
Since the number of reticulocytes varied from 100 to 500 per 1000 
erythrocytes, the variance according to the binomial law should vary 
from 90 to 250. This result led to a cooperative study [4] which in- 
cluded a control experiment with photographic technique. Our data 
also showed an extraordinarily strong subnormal variation. We hope 
that experiments now in progress will provide an explanation. 

Experimental design. In his work on biological standardization, 
Ipsen realized that mortality rates sometimes are relatively inefficient 
measures of response. He therefore recorded the survival time of the 
animals that died and developed a procedure for handling such data. 
By this means he increased by four-fold the efficiency of an experiment 
on the standardization of tetanus serum [3]. 

In experiments by Engbaek [1], groups of mice were given amounts 
of living culture of Bacillus Pfeiffer, type b, varying from 10** to 10°° 
bacteria. In this case neither the mortality rate nor the survival time 
gave a satisfactory reaction curve. The multiplication of the bacteria 
and their dispersal through the blood would seem to bear directly upon 
the development of the disease. Accordingly we studied the curve of 
bacteriemia in 18 mice by means of blood samples taken at frequent 
intervals. Our results indicated that the number of bacteria per cc. 
of blood should be determined routinely four hours after inoculation. 
This gave a steep reaction curve representing an increase in efficiency 
of at least 1000 per cent. 

The above two examples emphasize the importance of a careful pre- 
liminary investigation for selecting a suitable experimental object. 


All references for Recent Biometric Devel- 
opments in Denmark are listed on page 189. 
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RECENT WORK ON ‘‘INCOMPLETE BLOCK DESIGNS”’ 
IN INDIA 


R. C. 


The main principles of ‘‘experimental design’’ as it is understood 
at present were developed by R. A. Fisher and his associates at Roth- 
amsted in England during the years immediately following the first 
world war. The three basic principles of ‘‘experimental design’’ are 
‘‘randomization’’, ‘‘replication’’ and ‘‘loeal control’’. As Fisher has 
remarked, design and analysis are two aspects of the same logical whole, 
and randomization and replication are necessary in order that valid 
tests of significance based on the use of ‘‘t’’ and ‘‘z’’ tests should be 
possible. But to increase the efficiency of the experiment, ‘‘local con- 
trol’’ is necessary. The main principle is to divide the experimental 
material into relatively homogeneous parts called ‘‘blocks’’, and to so 
arrange the analysis that in making treatment comparisons ‘‘block 
effects’’ might be eliminated. Since any difference in the responses 
of the experimental units within the same block goes to swell the error 
(which serves as a measuring rod for the significance of the treatment 
differences) the experiment tends to become inefficient when the block 
size is large. We have, therefore, to find a way of reducing the block 
size. This has given rise to the device of ‘‘confounding’’ in the case 
of factorial designs, and the use of ‘‘incomplete block designs’’ in the 
ease of trials with a large number of varieties. I shall deal with some 
of the recent work done in India on the later class of designs. 

The incomplete block designs most commonly in use are the ‘‘lattice 
designs’’ and the ‘‘ balanced incomplete block designs’’. In the latter 
‘*y’’ treatments have to be arranged in ‘‘b’’ blocks of ‘‘k’’ experi- 
mental units each, so that every treatment is replicated ‘‘r’’ times, and 
each pair of treatments occurs together in A blocks, (the number A 
being independent of the pair of treatments we start with). It is 
readily seen that 


bk=vr, A(v-1) =r(k-1) 
and Fisher has shown that 


Given the parameters, v, b, r, k, A, satisfying the above conditions, 
the combinatorial problem of actually producing the required arrange- 
ment still remains. The practically useful cases are when r = 10, and 
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a number of these were listed in Fisher and Yates’ Tables. But many 
blanks remained. Indian workers have succeeded in filling up most 
of these blanks. Thus solutions have been found for the cases 


v=25, b=50, r= 8, k= 4, A=1 


v=41, b=82, r=10, k= 5, A=l 
y=28, b=63, r= 9, ke Awl 
v=25, b=25, r= 9, k= 9, A=3 
v=31, b=31, r=10, k=10, A=3 


v=21, 6=30, r=10, k= 7, A=3 
and it has been demonstrated that the following cases are impossible 


v=22, b=22, r= 7, k= 7, A=2 
v=16, b=21, r= 7, k= 5, A=2 
v=29, b=29, r= 8 k= 8 A=2 
v=21, b=28, r= 8, k= 6, A=2 


Some of these new solutions were included in the second edition of 
Fisher and Yates’ Tables and others are going to be included in the 
forthcoming edition. Only four cases with r=10 yet remain un- 
known. These are 


v=46, b=46, r=10, k=10, r=2 
v=36, b=45, r=10, k= 8 A=2 
v=46, b=69, r= 9, k= 6, A=l 


The total number of available balanced incomplete block designs is 
unfortunately very restricted. To a certain extent this deficiency is 
made good by Lattice designs in which, though every pair of treat- 
ments is not compared with the same accuracy, the analysis is simple 
and straightforward. In order to fit in with various experimental re- 
quirements which arise in practice, a much more flexible system of 
designs is required. Such a system can be obtained by the following 
considerations. It can be shown that any incomplete block experi- 
ment in which the block size k is constant, the effects can be estimated 
by solving the following system of normal equations 


t+ t+... 4+ Cinty=%, @=1,2,... 4, 


where v is the number of treatments and 9 is the ‘‘adjusted yield’’ 
corresponding to the 7th treatment, (i.e. 9; is the result obtained by 
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subtracting from the total yield of the 7th treatment, the average yield 
per plot of the blocks in which the treatment occurs). Also 


1 
cu=r(1-Z), 


where r; is the number of replications of the i-th treatment, and Ai; is 
the number of blocks in which the i-th and j-th treatments occur to- 
gether. In general, therefore, a system of v linear equations has to 
be solved, in order to estimate treatment effects. This is impracticable 
when v is large. We have, therefore, to impose combinatorial condi- 
tions, which should enable us to obtain an easy solution of the normal 
equations. One such set of conditions is the following: 


(1) Each block is of constant size and each treatment is replicated 
the same number of times. . 

(2) With respect to any treatment 6 the remaining ones can be 
divided into m groups of 7, M2, ... %m each, so that treatments of the 
+th group occur with the given treatment just A; times. (The numbers 
My, M2, ... Mm, Ar, Ac, .. . Am are constants of the designs, and are inde- 
pendent of the treatments with which we start). The treatments of the 
i-th group are said to be 7-associates of the original treatment. 

(3) The relation of association is symmetrical. Thus if the treat- 
ment 6 is an i-associate of , then @ is an 7-associate of 6. If pj,‘ is the 
number of treatments which are both j-associates of 6 and k associates 
of %, then the number of pj,‘ is independent of the pair of ?-associates 
@ and @ with which we start. 

Designs satisfying the above conditions have been called partially 
balanced incomplete block designs and include as special cases the 
‘‘balanced incomplete block designs’’ and the ‘‘lattice designs’’. The 
treatment differences are estimated with m different accuracies, and the 
solution of the normal equations can be made to depend on a solution 
of m simultaneous linear equations. Of course the practically useful 
cases are when m is small, say m = 3. 

Though a large number of illustrative cases have been worked out, 
a complete listing of the practically useful cases together with their 
combinatorial solutions remains to be undertaken. I am gratified to 
find that in the United States incomplete block designs, especially lat- 
tices, are being used on a large seale. and hope that before my de- 
parture from this country some of the new designs found by Indian 
workers may be actually tested in the field. 
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BIOMETRIC WORK IN AUSTRALIA 


Maurice BEuz 


In Australia the greater part of systematic biometric research is 
carried out by the Council for Scientific and Industrial Research 
through its various laboratories and institutions in the several states. 
This Council recruits its statisticians, in general, from University grad- 
uates who have had basic training in mathematics and, preferably, ad- 
vanced training in mathematical statistics. The Council passes these 
recruits through a preliminary period of practical training and then 
attaches them to the appropriate research department. 

University training in statistics at the uadergraduate level is pro- 
vided in most of the Australian universities. In Melbourne, for ex- 
ample, there are three courses available: (1) a pass course in theory of 
statistics, taken in the second year and following a fairly complete first 
year course of calculus, in which the standard distributions and tests 
are fully described and illustrated; (2) an honor course in mathe- 
matical statistics, taken either by mathematical graduates or by honor 
undergraduates who have completed a course in the theory of func- 
tions ; (3) a descriptive course in statistical methods for research work- 
ers, mainly graduates in bacteriology, biochemistry, pharmacology, 
physiology, ete., many of whom have had a basic training in algebra 
and the elements of calculus. This training is usually acquired at 
school, a special subject to meet the needs of such people being included 
among the matriculation group of subjects. All the above courses are 
increasing steadily in popularity. 


BIOMETRIC DEVELOPMENTS IN THE NETHERLANDS 
M. G. NEURDENBURG 


I regret that I am unable to review the situation of the statistical 
field in the Netherlands. The main purpose of my visit is to learn, 
instead, of the things you are doing in the United States. As a med- 
ical man, I would underline the necessity of thinking biologically. 

Statistics on infant welfare, hospitalization, ete., are basic for the 
organization of practical work such as of social security. Agricul- 
turalists, too, are working seriously on statistics in Holland. 
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DISCUSSION 


D. Nanda. I shall give a brief view of the plant-breeding work 
which is being done in India, much of it by Dr. V. G. Panse at the 
Institute of Plant Industry at Indore. 

Plant-breeding requires a study of the variability in the existing 
material. The fact that environmental variability is not of use to the 
plant-breeder and that only genetic variability leads to the improve- 
ment of a crop requires a separation of the genetic and non-genetic 
variation. The regression of the progeny means upon the parental 
values has been determined to give the proportion of the genetic to the 
total variation. This method is often used as a measure of the genetic 
variability and if the material is rich enough, the breeder initiates some 
selection. The selection program is carried out mostly with the help 
of the progeny row technique. 

If the material lacks genetic variability, the breeder hybridizes to 
create variability. The progenies from different crosses are quite 
heterogeneous due to the different combination of genes. The vari- 
ability so created is both genetic and non-genetic. If the material is 
grown in a randomized and replicated trial the variation can be sepa- 
rated into different parts. Such a study also reveals the effects of 
hybrid vigor and heterosis. The equations for this problem were given 
by Professor Fisher when he visited India in 1938. Further, several 
experiments are being performed to compare the yielding capacity of 
pure lines and crosses at various stages of propagation; i.e. F1, F2, Fs, 
ete. 

Another technique in common use is to select material on the basis 
of a discriminant function. This method was given by Fairfield 
Smith and is being developed further by Panse. The genetic advance 
achieved by selecting a certain percentage on the basis of discriminant 
function is subject to sampling error. Bartlett gave a method in 1939 
(Supplement to The Journal of the Royal Statistical Society) for 
determining the standard error of the genetic advance. I have tried 
to develop this theory further and have obtained the standard errors 
of the genetic advance and of the regression coefficients of the dis- 
criminant function model. 


N. Rashevsky. I wish to review briefly the discussions of this after- 
noon’s session. Our colleagues from the Old World have done amazing 
work in their respective countries and we can only bow in admiration 
to them for having achieved so much under such adverse conditions. 
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Let us look, however, at the scope of the papers presented and dis- 
cussed today. They are all strictly statistical in their nature. But 
an application of mathematics to any natural phenomenon goes far 
beyond the application of statistical method. The criterion of any sci- 
entific achievement is predictability. Frequently by the use of statis- 
tics we can predict rather precisely the future course of some events, 
just as by observing carefully the motion of a-golf ball we might pre- 
dict whether it will land in a hole or not. But you all will agree that 
there is a difference between this type of prediction and those exem- 
plified by Maxwell’s prediction of electromagnetic waves or by Ein- 
stein’s prediction of the transmutability of matter into energy. These 
latter predictions are made by first formulating some quantitative 
hypotheses concerning the intimate mechanisms underlying a set of 
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observable phenomena. By means of mathematical analysis conse- 
quences are drawn from this set of hypotheses and compared with 
available data. If the hypothesis is a good one, not only will its con- 
sequences agree with known facts, but they will lead to the prediction 
of new hitherto unknown phenomena. 

I wish to illustrate with a few figures the application of this method 
to biology, as it has been carried out in the last twelve years at the 
University of Chicago, although time does not permit me to explain 
the hypotheses. By making certain assumptions about the mechanism 
underlying cell division, mathematical equations describing the elonga- 
tion and constriction of freely dividing cells have been derived by Pro- 
fessor H. D. Landahl. To test those equations experiments were made 
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by Buchsbaum and Williamson here at Woods Hole. A typical case of 
agreement between the predicted course of the phenomenon and that 
observed is shown in Figure 1. I wish to emphasize that what you see 
here is not an empirical fitting by statistical methods of a curve to an 
observed set of points. The analytical expression for these curves was 
obtained long before the experiment was planned. 

Professor Landahl and others have also contributed to the theo- 
retical study of the functions of the central nervous system. By mak- 
ing some rather complex assumptions about the nature of interaction 
of the neurons in the brain, conclusions can be drawn, for instance, as 
to the dependence of the reaction time to a stimulus on the interval 
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between that stimulus and a ‘‘preparatory’’ stimulus given as a warn- 
ing in advance. The comparison of the theory and the experiment, 
taken also from Landahl’s publication, is shown in Figure 2. 

Figure 3 shows the agreement between theory and experiment on 
the relation between the intensity of a stimulus and the just noticeable 
difference between two stimuli which still can be discriminated. These 
are typical curves taken from a series of papers by Householder. 

I have no time to discuss a number of other important contributions 
in the field of mathematical biophysics such as Landahl’s studies on 
learning, Weinberg’s work on excitation, Williamson’s work on elec- 
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trical properties of cells, Bloch’s work on permeability, Sacher’s work 
on periodicities in neural circuits, Morales’ work on gas exchanges, and 
many others. I also cannot mention the innumerable questions to 
which we still have no answers. However small our knowledge and 
however great our ignorance, I repeat that only through the har- 
monious simultaneous development of all possible approaches to mathe- 
matical biology shall future progress be made. 
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Saturday morning, September 6 


This session, under the chairmanship of Drs. Teissier and Dieule- 
fait, dealt with international cooperation in Biometrics. The first 
speaker was Dr. D. W. Bronk of the National Research Council. Dr. 
Bronk emphasized the importance of free communication between the 
scientists for all nations in establishing a stable peace, and discussed 
the various forms of international organization which have been 
adopted in the past. 

Dr. Dieulefait then spoke of the reorganization of the International 
Statistical Institute. Under its new statutes it would include not only 
governmental and economic statistics but also mathematical statistics 
and its applications. He suggested that the society which might be 
formed this morning should consider affiliation with 1.8.1. 

Professor R. A. Fisher then submitted the report of the Committee 
on International Organization, which recommended the formation of 
an international membership society. A draft constitution was dis- 
tributed in mimeographed form and considered article by article under 
the chairmanship of Dr. Maurice Belz. 

The inclusion of the word ‘‘international’’ in the name of the 
Society was debated first. It was finally agreed that the name of the 
organization shall be the ‘‘ Biometric Society’’ and that a short descrip- 
tion clarifying the objectives of the Society shall appear on its letter- 
heads. After further small modifications, the draft constitution was 
adopted. 

It was then resolved that all persons present at the Conference, or 
invited but unable to attend, shall be Charter Members of the Biometric 
Society if they so desire. It was further resolved that the Committee 
on International Organization and such other Members as it may 
designate shall constitute the first Council of the Society. 

The Conference, then sitting as the Biometric Society, unanimously 
adopted the constitution as amended during the morning session. A 
copy of this constitution follows. 
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THE BIOMETRIC SOCIETY—CONSTITUTION 


1. Scope of the Society 


The Biometric Society is an international society for the advance- 
ment of quantitative biological science through the development of 
quantitative theories and the application, development and dissemi- 
nation of effective mathematical and statistical techniques. To this 
end the society welcomes to membership biologists, mathematicians, 
statisticians and others interested in applying similar techniques. 


2. Members 


The Society shall have one class of members. To become a member, 
a person must be proposed by two members of the Society and ap- 
proved by the Council. The Council may delegate this authority. The 
members represent the highest authority of the Society. The Council 
shall consult them on any vital questions that may effect the policy of 
the Society as a whole, obtaining their decision by a mail vote. 


3. Officers 


The general officers of the Society shall be the President, the Secre- 
tary and the Treasurer. The regional officers shall be the Vice-Presi- 
dents, Regional Secretaries and Regional Treasurers. The officers shall 
be elected by the Council for one year. 

The President shall act as chairman of the Council. 

Each Vice-President shall represent a specified region. The regions 
shall be determined by the Council from time to time. 

The Treasurer shall present financial statements to the Council and 
shall bring condensed statements to the attention of the members. The 
offices of Secretary and Treasurer may be combined. 

No President or Vice-President shall serve more than two consecu- 
tive years. 

The Council shall elect Regional Committees to serve under the 
chairmanship of each Vice-President and may elect Regional Secre- 
taries and Regional Treasurers where that seems advisable. 

The general and regional officers, acting as a nominating committee, 
shall submit to the Council a list of candidates for general and regional 
officers but members of the Council may vote for names not on this list. 
Regional officers and committeemen may be nominated directly by the 
members at regional business meetings. 
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4. The Council 


The President, the Vice-Presidents, the Secretary, and the Trea- 
surer shall be ex-officio members of the Council. There shall not be 
less than twelve nor more than twenty ordinary members of the Coun- 
cil, who shall be elected with a view to representing the various geo- 
graphic areas and fields of activity in which the Society has members. 
Ordinary members of the Council shall be elected for three years by a 
mail vote of the members, the terms of approximately one-third of them 
terminating each year. 

The general and regional officers of the Society, acting as a nomi- 
nating committee, shall submit to the members a list of names of twice 
the number of candidates necessary to fill the ordinary vacancies in the 
Council, but members may vote for names not on this list. No ordi- 
nary member of the Council shall serve more than two consecutive 
terms of three years each. 

In matters of purely regional importance, the appropriate Regional 
Committee may act for the Council, but the Council may reverse its 
action. 

5. Activities 


Any activities which fall within the sphere of interest of the Society 
may be authorized by the Council, such as international and local sci- 
entific meetings and the issuance of publications reporting the activities 
of the Society or containing other matters of biometric interest. 

The Society may affiliate itself with international organizations. 
With Council approval a Region may affiliate itself with regional or 
national organizations. 


6. Financial Organizations 
The dues for members shall be fixed by the Council. Requests and 
gifts may be received. 
7. Amendment of the Constitution 


Amendments to the Constitution must be approved by the Council, 
or in writing by at least five per cent of the members, and before be- 
coming effective they must be ratified by a two-thirds majority of those 
voting in a mail vote taken among all members. 
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Saturday afterncon, September 6 


The Chairman, Leslie F. Nims, called upon Professor Wilson, who 
presented the following resolutions. 


“The members of the First International Biometric Conference held 
at Woods Hole, September 5-6, 1947, desire to express their profound 
thanks to the Marine Biological Laboratory, to the Rockefeller Founda- 
tion, to the Connecticut Agricultural Experiment Station, to the National 
Research Council and to the Officers and Organizing Committee of the 
Conference for providing facilities and making arrangements which have 
resulted in an initial Conference of great present success and high hope 
for the future.” 


It was voted unanimously that the Secretary be directed to transmit 
this resolution to all parties concerned. The Chairman then intro- 
duced Professor Teissier, who gave his paper* on ‘‘La Relation d’Allo- 
metrie, sa Signification Statistique et sa Logique’’. Dr. J. Monod 
opened the discussion by summarizing Dr. Teissier’s address in English. 

The First International Biometric Conference was concluded with 
the following remarks by Professor Teissier. 


ALLOCUTION DE CLOTURE DE LA PREMIERE CONFER- 
ENCE INTERNATIONALE DE BIOMETRIE 


GEORGES TEISSIER 


Avant que nous nous séparions, je voudrais, mes chers collégues, 
vous remercier encore de l’honneur que vous m’avez fait en me 
chargeant de présider vos débats, et vous renouveler mes excuses pour 
le retard bien involontaire qui m’a privé du plaisir d’ouvrir votre con- 
férence et d’entendre ]’exposé du Professeur Fisher. 

Je voudrais aussi remercier en votre nom 4 tous ceux qui ont eu 
l’initiative de cette réunion et qui ont réussi 4 l’organiser si heureu- 
sement. C’est tout particuliérement au Dr. Bliss que va notre recon- 
naissance. Nous savions, dés longtemps, ce que lui doit la Biométrie 
et je ne suis pas prés d’oublier combien ses méthodes de calcul m’ont 
été personnellement utiles en de certaines circonstances. Nous sommes 
maintenant quelques dizaines de biométriciens 4 savoir, en outre, qu’il 
est aussi un organisateur habile et le plus sympathique des hétes. Qu’il 
recoive ici publiquement 1]’assurance de notre gratitude 4 tous. 


_ * This paper and the ensuing discussion will be published in the March issue of 
Biometrics. 
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Puisque les fonctions que vous avez bien voulu m’attribuer me 
donnent le droit de parler en dernier, je me permettrai de vous 
présenter pour conclure quelques réflexions sur notre nouvelle asso- 
ciation. 

Au point ot en sont arrivées, dans leur développement, la Statis- 
tique et la Biologie, je crois qu-une tentative telle que celle que nous 
réunit ici était devenue nécessaire, mais je crois aussi que l’entreprise, 
dans laquelle nous nous lancons d’un méme élan, est difficile. Elle ne 
réussira pleinement que si, dés le début, nous avons une conscience 
claire de la nature exacte de la tache qui nous attend. 

Notre Société ne doit pas étre une nouvelle association de statisti- 
ciens et pas davantage une nouvelle association de biologistes. Nous 
appartenons les uns et les autres a de telles sociétés et nous continue- 
rons, sans doute, a participer de notre mieux a leurs travaux. Dans la 
Société de Biométrie, il nous faudra travailler ensemble a faire progres- 
ser la Biologie par la Statistique. Pour cela, il nous faudra d’abord 
bien nous comprendre. Certains d’entre nous ont surtout une forma- 
tion mathématique et connaissant assez mal les problémes biologiques ; 
d’autres surtout biologistes n’ont pas poussé bien loin leurs études 
mathématiques. Bien peu nombreux sont ceux qui ont eu le bonheur 
de bénéficier, dés leur jeunesse, de cette double culture mathématique 
et biologique qui apparait cependant comme essentielle 4 la pleine 
compréhension de la science que nous cultivons. 

Il nous faudra donc—et ce devra étre notre premier objectif—nous 
instruire mutuellement, le statisticien apprenant au contact du biolo- 
giste quelles sont les préoccupations dominantes de celui-ci, le biologiste 
apprenant de son collégue statisticien |’art difficile de poser correcte- 
ment un probléme quantitatif. Il nous faudra ensuite recruter, parmi 
les jeunes, des futurs biométriciens que nous devrons former, mieux 
que nous ne |’avons nous-méme été, 4 cette passionnante mais sévére 
discipline. 

Le fait que nous soyons réunis aujourd-hui, si nombreux et venus 
de tant d’horizons divers, prouve que nous avons compris combien était 
utile au progrés de la science, cette confrontation d’hommes d’origine 
et de formation dissemblables, mais épris du méme désir de connais- 
sance désintéressée. Nous pouvons bien augurer de notre avenir 
commun, si nous savons intéresser un plus grand nombre de biologistes 
a nos efforts et si nous savons leur montrer quel merveilleux secours 
apporte a ]’étude des phénoménes de la vie, ]’instrument mathématique 
sagement manié. 
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La double organisation de notre association, 4 le fois internationale 


et régionale, nous aidera 4 assurer ce recrutement et c’est avec con- 
fiance que nous pouvons lancer a travers le monde la nouvelle Société 
de Biométrie. 


[2] 


[3] 


[4] 


[5] 


[6] 


[7] 


[8] 


(9] 


[10] 
[11] 


[12] 
[13] 
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THE BIOMETRIC SOCIETY 


REPORT OF COUNCIL MEETINGS 


First meeting of the Council, Woods Hole, September 6, 1947, with 
the following present: Belz, Bliss, Bose, Fisher, Hopkins, Linder, 
Teissier, Rasch and Tukey with Monod as interpreter and guest. As 
the first order of business, Professors Cox, Haldane and Wilson were 
elected members of the Council. Since she was still in Woods Hole, 
Miss Cox joined the Council for the remainder of the meeting. As its. 
first officers the Council elected R. A. Fisher, President ; J. W. Hopkins, 
Treasurer; and C. I. Bliss, Secretary. 

The publication of the Proceedings of the Woods Hole Conference 
and the policy of the Society in respect to a journal were discussed. 
As editor of Biometrics, Miss Cox offered to publish the Proceedings 
in the December issue (1947) and the two main papers by Professors 
Fisher and Teissier with the ensuing discussions in the March 1948 
issue. After considering an alternative proposal of The American 
Naturalist, the Council agreed to accept Miss Cox’s offer. It was 
agreed that the papers from the Proceedings should be published in 
the language in which they were delivered. 

Members of the Council favored an arrangement with the American 
Statistical Association, by which Biometrics could be used as the offi- 
cial journal of the Biometric Society during 1948. It was agreed that 
dues should include a subscription to Biometrics, except that the mem- 
bers already receiving the journal through the Biometrics Section 
would receive a rebate. Before adjourning, the Council agreed to 
meet next in Washington. 

Second meeting of the Council, Washington, D. C., September 15, 
1947, with the following present: Belz, Bliss, Bose, Cox, Dieulefait, 
Fisher, Linder, Neurdenburg, Rasch and Tukey. The following were 
elected to the Council: Buzzati-Traverso, Cole, Demerec, Goulden, 
Johnson, Needham and Simpson. 

It was decided to set up a skeleton regional organization but allow 
members the option of belonging to a given region or to be members at 
large. Four regions were set up with officers: British Region, Indian 
Region, Western American Region (including the American Pacific 
coast, western Canada and Mexico) and Eastern American Region (in- 
cluding eastern Canada). The following five regions were to be acti- 
vated as soon as it seemed desirable: Scandinavian Region (Denmark, 
Norway, Sweden and Finland), Benelux Region (Belgium, Netherlands 
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and Luxemburg), Australian Region, French Region and Russian 
Region. Each member of the Council agreed to handle publicity for 
the new Society in his region. 

As an incentive to membership it was decided to consider everyone 
joining before February 1, 1948, as a charter member of the Society. 
The Secretary was given instructions by the Council, which would 
facilitate the election of additional charter members. 

In accord with Article 5 of the Constitution, the Council approved 
affiliation of the Society with the International Council of Scientific 
Unions if such could be arranged. President Fisher was authorized 
to review the matter with the General Secretary of the 1.C.S.U. The 
Council also approved affiliation with the International Statistical In- 
stitute, when its statutes make this possible. The Eastern American 
Region was encouraged to affiliate with the National Research Council’s 
Division of Biology and Agriculture, with the American Association 
for the Advancement of Science and with the American Statistical 
Association. 

In a discussion of finances the Secretary reported that a balance of 
possibly $300 from the Rockefeller grant in support of the Woods 
Hole Conference would be available for printing and distributing the 
Proceedings of the meetings of September 5th and 6th but not for the 
use of the new Society. The Board of Directors of the American 
Statistical Association has offered to spend up to $100 for printing 
extra copies of the December issue of Biometrics for distribution to 
prospective members of the Biometric Society who would become sub- 
seribers. In consideration of the probable cost of Biometrics to the 
Society in bulk subscriptions, the Council voted dues of $4.00 for 1948 
with a rebate of $2.00 to those who are also members of the Biometrics 
Section. The policy was approved of having dues collected by the 
regional secretary-treasurers who would be authorized to retain $1.00 
or its equivalent for regional use. In view of the limitations upon 
international exchange, dues are to be payable in local currencies 
except in an amount sufficient tc cover a subscription to Biometrics, 
which must be received in dollars before the journal can be mailed. 

The letterhead of the Society is required to carry a descriptive 
clause beneath its name. After some discussion, the Council decided 
on the wording ‘‘ An International Society Devoted to the Mathemati- 
cal and Statistical Aspects of Biology.’’ 

C. I. Buss 
Secretary 
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ABSTRACTS 


(42) 

UPHOLT, WILLIAM M. AND SCUDDER, H. I. (U. S. Public 
Health Service.) A Problem in Sampling Mobile Insect Populations. 

Some insect population distributions are so extremely ‘‘con- 
tagious’’ that a reasonable number of random samples does not 
provide a satisfactory estimate of the mean. In certain cases such 
as house flies the insect itself is so highly mobile that no satisfactory 
method of obtaining a truly random sample has been devised. To 
meet this situation one of the authors has suggested a sampling 
method which apparently provides an estimate of the magnitude of 
the greatest concentrations of flies. These estimates appear to be 
reproducible and to have a relatively small, though undetermined, 
variance. The problem presented is the statistical validity of, and 
the selection of proper statistical techniques to use with, such 
samples which appear to be empirically justified from the stand- 
point of field biologists. 


(43) 

COOK, 8. F. (U. C. L. A.). Survivorship in Aboriginal Populations. 

Survivorship in truly aboriginal populations can be determined 
by only two methods: (1) censuses taken by members of an invad- 
ing culture prior to any disturbance caused by such invasion and 
(2) determinations of the age at death of skeletal remains. It has 
been possible to assemble data concerning several such populations, 
including Tasmanians, East Africans and a number of North 
American aboriginal groups. The mean age at death varies widely. 
Thus the Pecos Indians (Hooton) and those of the Valley of Mexico 
appear to have survived to a relatively old age whereas the Tasma- 
nians (Todd) and the California Indian tribes tended to die at a 
very early age. The conclusion is reached that survivorship is not 
an inherent or genetic characteristic but is dependent purely upon 
environmental conditions and habits of life. Moreover the data 
studied appear to indicate that in peoples wholly untouched by out- 
side influences infant mortality is less than has been commonly sup- 
posed. 
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QUERIES 


(54) 
QUERY: This is purely a computational matter, but one which has 
bothered me frequently, and upon which I should like to get ad- 
ditional information if there is any available. In many calcula- 
tions I multiply a series of quanties 2, 22, 23, . . . tn by a con- 
stant factor k, rounding the individual products to the nearest 
integer. As a check I compare Skx (rounded) with k(Sz). 

Now the difficulty is that I have no idea of the expected distribu- 
tion of the difference Ska—k(Sx), beyond the obvious facts that it 
has minimum 0 and maximum »/2 (disregarding signs). Through 
long experience I have arrived at the conclusion that where n is 
around 20-50, a value of Ska—k(Sz) of over 2 usually indicates an 
error in one or more of the individual products. I should like to 
treat the matter more exactly if possible, and wondered if you were 


aware of anything on the subject in the various works for com- 
puters, or elsewhere. 


ANSWER: Your problem is that of the distribution of means (or 
totals) of samples from a rectangular population. This  distribu- 
tion is diseussed in The Advanced Theory of Statistics, Volume 1, 
by M. G. Kendall on pages 240-242. 
If we denote the total of sample values from this distribution by 
z, the distribution of z has the mean zero and the variance »/12. 
Since the distribution of z tends to normality with increasing n, for 
moderately large values of » we may ascribe to a total a standard 
error of \/n/12. The difference between S(kz) and k(Szx) should 
exceed twice this standard error in only about 5 percent of cases. 
If n is, for example 50 the standard error of the total will be about 
2 (that is, \/50/12). If the discrepancy is greater than 4, there is 
good reason to suppose that some error of calculation has been made. 
Exact levels of significance have not so far as I know been evaluated, 
but the above rule should be reasonably accurate. 
Oscak KEMPTHORNE 
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(55) 
QUERY: This is to request information as to the applicability of 
the analysis of variance procedures and the F-test to the studying of 
the significances between classifications of variance rather than 
means. That is, in a table of variances (8?) as follows, could the 
variance of the variances by rows, by columns, and a residual (in- 
teraction) be computed following the usual procedure and then by 
the F-test could the significance of the differences in the mean vari- 
ances for the columns or rows be determined ? 


S*,, S*,, S*,, 
8? 12 ‘ S722 S? 32 


S23 


ANSWER: Bartlett and Kendall discuss this problem in ‘‘The 
Statistical Analysis of Variance—Heterogeneity and the Logarith- 
mic Transformation’’, Supplement to the Journal of the Royal 
Statistical Society, Vol. 8, pages 128-138, 1946. If you have inde- 
pendent samples of approximately the same size, if they are drawn 
from normally distributed populations, and if they have ten or more 
degrees of freedom, you can transform to the logarithms of the 
variances and analyze the variance of the resulting logarithms in 
the usual manner. Evexu if the degrees of freedom in the samples 
are as few as five, the method may be used tentatively. 

If you need information as to why transformation is necessary, 
you will find excellent discussions of analysis of variance by Eisen- 
hart, Cochran, and Bartlett in the March issue of Biometrics. 

Other articles describing approximate tests of homogeneity of 
variances are as follows: Stevens in the Journal of Genetics, Vol. 33, 
page 398; Cochran in the Supplement to the Journal of the Royal 
Statistical Society, Vol. 4, page 102; and Bartlett in the same Jour- 
nal, Vol. 4, page 137. 


Water T. FEDERER 
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NEWS AND NOTES 


G. E. DICKERSON, formerly with the Regional Swine Breeding 
Laboratory at Ames, Iowa, has a new address—Animal Husbandry 
Department, University of Missouri, Columbia. He is in charge of 
research and teaching in animal breeding. Animal breeding there 
includes animal genetics and physiology, particularly physiological 
reproduction. He hopes to be conducting genetic studies of farm 
and laboratory animals ... JAMES H. BYWATERS, Virginia Agri- 
cultural Experiment Station, Blacksburg, expects to plan and con- 
duct poultry research projects aimed toward development and guid- 
ance of the poultry industry in Virginia. Welcome to the South! 
Mr. Bywaters comes to Virginia from the U. S. Department of Agri- 
culture, Poultry Research Laboratory at East Lansing, Michigan 
...d. F. KaGy, Agricultural Research Laboratory, The Dow 
Chemical Company, Seal Beach, California, writes, ‘‘It is true that 
most of our populations are skewed and if we were to obtain ac- 
curate evaluations of our experimental errors and treatment dif- 
ferences, transformations would have to be used. We are greatly 
in need of suggestions and help concerning the problem of experi- 
mental design. This problem is the single greatest difficulty with 
which many experienced research men are confronted. The tend- 
ency is to make experiments too complicated.’’ Quite true. We 
hope to present a series of articles dealing with the designing of 
experiments in the near future. There have been numerous re- 
quests for illustrations from actual experiments . .. ANTONIO E. 
MARINO, Sarmiento 246, Buenos Aires, Argentina, writes, ‘‘I remem- 
ber the nice days I spent at Raleigh in 1941. At present I am work- 
ing with the Argentine Seed Branch of Cargill Incorporated as a 
plant breeder, specially for developing corn hybrids for this coun- 
try.’’ ... FERNANDO VILLAMIL-GaRcIA, Estacion Agricola Experi- 
mental, Palmira-Valle, Colombia, South America, has had a few 
‘solos’? into the statistical world even though he struggled a little 
when taking Mathematics 441 and 442 with Professor Snedecor! 

. F. W. PRESTON, Preston Laboratories, Butler, Pennsylvania, 
states, ‘‘although we are not concerned very much with statistics as 
applied to biology, we find that biologists tend to develop branches 
of statistics that are useful to us in glass technology. Here, unlike 
most branches of physics, we run into wide fluctuations of our 
measured quantities, and the simpler forms of statistics are not 
always adequate’... c. BIRCH, Australian visiting research 
worker at the Bureau of Animal Population, Oxford University, 
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thinks that Biometrics can continue to serve its most useful function 
by summarizing some of the new developments in statistics... . 
H. P. DONALD, Institute of Animal Genetics, University of Edin- 
burgh, requests that Biometrics contain papers which deal with 
practical problems in a way intelligible to people with little mathe- 
matical skill. We would like to meet these objectives, but where are 
enough biological statisticians who can take time to write for 
Biometrics! . . . In discussing some of these points ERIC C. WOOD, 
Virol Limited, Ealing, England, recently made a profound state- 
ment in a letter to the editor. ‘‘An experimental design which is 
both sound in theory and convenient in practice is usually evolved 
only by two people, a statistician and a practising assayist in col- 
laboration.’’ . . . GEORGE W. HARVEY has been appointed director of 
the Statistical Division of the newly organized National Blood Pro- 
gram of the American Red Cross. He will institute modern pro- 
cedures for the collection and analysis of relevant data, pertaining 
to blood collection and the fractionation of blood. This new pro- 
gram will provide blood and blood derivatives throughout the 
nation to help save lives, prevent needless suffering, and to further 
research. . . . ERRETT C. ALBRITTON, Professor of Physiology, The 
George Washington University School of Medicine, Washington, 
D. C., feels that some periodical might give space to the analyses of 
problems taken from current literature. He states, ‘‘The medical 
investigator has nowhere available, as a model for his own thinking, 
such explicit logical analyses of the research of others.’’ A statis- 
tical review of medical research sounds difficult! . .. F. BE. ALLISON, 
Senior Chemist, Bureau of Plant Industry, Beltsville, Maryland, 
has expressed an idea we have heard from others. ‘‘I might add 
that usually a research worker must be his own statistician because 
there is no one near with whom he can consult. Often he is greatly 
puzzled as to what procedure to follow. A book that gives a num- 
ber of typical experimental designs, suitable for various types of 
problems, and sets forth in detezil the methods of calculation to be 
used, together with the logical conclusions to be drawn from the 
data after analysis, would, it seems to me, be of very great value. 
The information is now available, of course, but not readily, and 
certainly not in a single volume. Such a book on experimental 
designs and analytical procedures should not be limited to analysis 
of variance methods but should include other common procedures of 
use in obtaining or analyzing data. To be of most use the book 
should be as simple and non-mathematical as possible.’”?” ... F. x. 
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LAUBSCHER, senior research officer, College of Agriculture and Ex- 
periment Station, Potchefstroom, Union of South Africa, is in this 
country on a commission to investigate fiber production and to visit 
experiment stations where corn breeding is being done. At Raleigh, 
he expressed interest in the statistical approaches to breeding tech- 
niques of peanuts, soybeans and small grains. . . . The Statistical 
Section of the Tennessee Public Health Association met May 6 (its 
first anniversary) in the Andrew Jackson Hotel, Nashville. They 
held both a morning and afternoon session with the following program: 
Morning Session—9:00 A.M.—1. Program of the Medical Research 
Statistics Division of the Veterans Administration—PAUL M. DENSEN, 
Chief of Division of Medical Research Statistics, Veterans Adminis- 
tration. 2. Statistical Needs of Social Workers—LORA LEE PEDERSON, 
Director of Nashville School of Social Work. 3. Fisher’s Exact Sig- 
nificance Test for the Fourfold table—mARGARET MARTIN, Assistant 
Professor of Preventive Medicine and Public Health, Vanderbilt 
University School of Medicine. *%4. Statistical Training Program of 
Tennessee Department of Public Health—sara LoU HATCHER, Statis- 
tician, Tennessee Department of Public Health. 5. Completeness of 
Data on Birth Certificates—p. c. PETERSON, Director of Division of 
Vital Statistics, Tennessee Department of Public Health. Afternoon 
Session—2 :00 P.M.—1. Public Health Statistical Service in Mississippi 
—MARGARET E. RICE, Supervisor of Public Health Statistics, Mississippi 
State Board of Health. *2. Source of Report of Notifiable Diseases— 
MRS. THOMAS PARRISH, Statistical Aide, Tennessee Department of 
Public Health. *3. Study of Reported Cases of Tuberculosis—LoLa 
MAI UPCHURCH, Statisician, Tennessee Department of Public Health. 
4. Observations Regarding Health Work in Chile—rUTH R. PUFFER, 
Director of Statistical Service, Tennessee Department of Public Health. 
Paul Densen presided. Mimeographed copies of the three papers with 
an asterisk are available upon request from the Tennessee State Depart- 
ment of Public Health. Thirty-four persons attended. The following 
officers were elected for the coming year: Chairman, ANN DILLON, 
Tennessee Department of Public Health ; Vice Chairman, HELEN CHAND- 
LER, Tennessee Valley Authority, Wilson Dam, Alabama; and Sec- 
retary, MARGARET MARTIN. 
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Power of test, 14, 22, 25 

Probit, 41, 47, 62 

Pure variance, 54, 77 

Quantitative analysis, 61 

Quartile, 56 

Randomization, 33 

Randomized blocks, 27, 40 

Random variable, 9 

Ranks, 49, 119, 138 

Rate, average, 56 

Reaction time, 172, 182 

Regression, 47, 66, 91, 103, 105, 137 

Rejection of data, 2 

Research in statintioal methods, 24 

Sampling, 19, 101, 163, 193 
systematic, 102 

Sampling error, 92 

Scores, 49 

Sensitivity of test, 22, 28 

Serology, 173 

Sex ratio, 99 

Skewness, 25, 

Source of variation, 92 

Spurious correlation, 65 

Square-root 41, 91 

Standard error, 100 

Statistical model, 163. 

T test, 24, 25, 28, 97 

Teaching of statistics, 1, 111, 179 

Test of normality, 26 

Test of significance, 3, 14, 24, 119, 
137, 139 

Transformation, 26, 39, 91, 98, 174, 
194 

pure, 54 

Weights, do, "4, 163 

X-ray dosage, 60 

Yates’ correction, 125 
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STALE OF THE OWNERSHIP MANAGEMENT, ETC., RE- 
QUIKE)D BY THE ACTS OF CONGRESS OF AUGUST 24, 191., AND MARCH 83, 1933, 
OF tries, published quarterly at \\ashington, D. C., for 12 months ending October, 
1947. 

Wasi .ngton, D. C., ss. Before me, A Notary Public in ind fcr the State and 
coun:, sforesaid, personally appeared Lester S. Kellogg, who, | aving been duly sworn 
accordin,s to law, deposes and says that he is the managing editor 0. Biometrics and 
that the following is, to the best of his knowledge and belief, a true -tatement of the 
ownerrh'9, management (and if a daily paper, the circulation), ete., of the aforesaid 
publication for the date shown in the .bove caption, required by the A‘t of August 24, 
1912, «s a-rended by the Act of Marc!: 3, 1923, embodied in section 5:7, Postal Laws 
and Regrl: tions, printed on the rever:e of this form, to wit: 

1. Tha’ the names and addresses of the publisher, editor, inana;:ing editor, and 
business m naagers sre: Publisher, A ierican Statistical Associa-ion, 1603 K Street, 
N. W., Was iington 6, D. C.; Editor, Gertrude M. Cox, Institute of Statistics. Raleich, 
N. C.:; Manising Editor, Lester 3. Kellogg, 1603 K Street. N. W., Washington 6, D. C.; 
Business M:nagers, None. 

2. That the owner is: (If owned by a corporation, its name and address must be 
stated and also immediately thereunder the names and addresses of stockholders owning 
or holding ove per cent or more of to‘al amount of stock. If not owned by a corpora- 
tion, the names and addresses of the individual owners must be given. If owned by a 
firm, company, or other unincorporated concern, its name and address, as well as those 
of each individual member, must be siven.) American Statistical Association, 1603 K 
Street, N. W., Washington 6, D. C. 

3. That the known bondholders, mortgagees, and other security holders owning or 
holding 1 per cent or more of total amount of bonds, mortgages, or other securities are: 
(If there are none, so state) None. 

4. That the two paragraphs next above giving the names of the owners, stock- 
holders, and security holders, if any, contain not only the list of stockholders and 
security holders as they appear upon the books of the company but also, in cases where 
the stockholder or security holder eppears upon the books of the compary as trustee or 
in any other fiduciary relation, the name of the person or corporation for whom such 
trustee is acting. is given: also thut the said two paragraphs contain statements embrac- 
ing affiunt’s full knowledge und belief as to the circumstances and conditions under 
which stockholders and security holders who do not appear upon the books of the com- 
pany as trustees, hold stock and securities in a capacity other than that of a bona fide 
owner; and this affiant has no reason to believe that any other person, association, or 
corporation has any interest direct or indirect in the said stock, bonds, or other securi- 
ties than as so stated by him. 

5. That the averaze number of copies of each issue of this publication sold or dis- 
tributed. through the mails or otherwise. to paid subscribers during the twelve months 
preceding the date shown above is (...... ). (This information is required from daily 
publications only.) 

Lester S. KELLocG, Mg. Ed. 

Swern to and subscribed beicre me this 2nd day of October, 1947. 

commission expires Jul, 31, 1249.) 


BrverLy RutH Rison, Notary Public 


» 
a 
le 
f 
: 


